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The nonlinear dynamics of axisymmetric, as well as helical, frozen-in vortex structures is investi- 
gated by the Hamiltonian method in the framework of ideal incompressible electron magnetohydro- 
dynamics. For description of current-sheet formation from a smooth initial magnetic field, local and 
nonlocal nonlinear approximations are introduced and partially analyzed that are generalizations 
of the previously known exactly solvable local model neglecting electron inertia. Finally, estima- 
tions are made that predict finite-time singularity formation for a class of hydrodynamic models 
intermediate between that local model and the Eulerian hydrodynamics. 
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Introduction. — It is a well known fact that current 
sheets play exclusively important role in plasma dynam- 
ics (see, e.g., |2| ^, ^ and references therein). However, 
analytical study of current sheets formation and their 
dissipative dynamics is a very difficult problem in the 
framework of usuaUy used nonlinear (and also nonlocal in 
the incompressible limit) equations of motion of plasmas. 
That concerns the usual magnetohydrodynamics (MHD), 
the electron magnetohydrodynamics (EMHD), as well as 
the multi-fluid models of plasmas. So, up today we do 
not have a mathematically clear answer on the question 
whether the current density will become singular in a fi- 
nite time or its growth can be only exponential in these 
systems. Numerical simulations remain to be the main 
tool for obtaining quantitative results ^, ||, ^ ||, |[ . 
Therefore, an important role for theoretical understand- 
ing of current sheets dynamics can be played by local 
nonlinear approximations, that sometimes have exact so- 
lutions describing formation of singularities. An example 
of such relatively simple, approximate differential equa- 
tion for the magnetic field B(r, i) is (see, e.g., [|| for 
derivation and explanation) 



Bt = curl 

Aire 
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This equation describes the motion of magnetic struc- 
tures in EMHD on length scales much larger than the 
inertial electron skin-depth, while the main part of the 
energy is concentrated in the magnetic field, with the 
kinetic energy of the electron fluid motion being much 
smaller. The equation (|^) has been extensively exploited, 
for instance, to study fast penetration of magnetic field 
into plasmas due to the Hall effect |l^, as well as 
rapid dissipation of magnetic fields in laboratory and as- 
trophysical conditions . The interest to this equation 
is explained, in particular, by the fact that axisymmetric 



configurations with B | 
able (see |, g @ 0) 



have been found exactly solv- 
. In this geometry, the equation 
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of motion is reduced to the well known one-dimensional 
Hopf equation, that should be solved independently for 
each value of the radial coordinate. The mechanism of 
singularity formation in these solutions is connected sim- 
ply with breaking in a finite time of the magnetic field 
profile. The magnetic field itself docs not become infinite, 
but its curl tends to the infinity at some point of the ax- 
ial cross-section. Inclusion of dissipative terms into the 
equation stops the breaking, but instead of multi-valued 
profile, a shock forms, the length of which increasing with 
time. The shock is a cross-section of a current sheet. 

The main purpose of present work is to extend the 
analysis of such axisymmetric flows by consideration ad- 
ditional nonlinear effects caused by electron inertia. They 
either play role of small corrections for long-scale flows, 
or, when shock becomes narrow, change drastically the 
dynamical behavior by smoothing the transport veloc- 
ity field. This situation is quite different in comparison 
with the self-similar EMHD solutions discussed in Ref. 
p^ . Also, the flows with other geometrical symmetry are 
considered below in the approximation (^), when all the 
frozen-in magnetic lines have helical shapes with a same 
spatial period along z-direction. In this case different 
level contours of the axial component of the magnetic 
field rotate in a perpendicular plane with different an- 
gular velocities, thus producing the shock. Finally, we 
use the developed technique to predict formation of fi- 
nite time singularities of the shock type in a class of hy- 
drodynamic systems that are in some sense intermediate 
between the model (0) and the usual Eulerian hydro- 
dynamics. That is a contribution towards general un- 
derstanding of possible mechanisms for singularities in 
hydrodynamics, that is a long-standing theoretical prob- 
lem. 

Incompressible two-fluid model. — Before the main 
consideration, it is useful to recall the place of EMHD 
among different hydrodynamical plasma models If 
there are only two kinds of particles in the plasma — neg- 
atively charged electrons with the mass m and positively 
charged ions with the mass M , then the most general is 
the two-fluid model, which contains MHD, EMHD, and 
Hall MHD as special cases. Let the equilibrium concen- 
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tration of particles of each sort be equal to n. If the tem- 
perature of the system is sufficiently large, nT ^ B^, 
then for slow vortical flows one can neglect deviations 
of the concentrations from n (the quasi-neutrality con- 
dition), and believe the velocity fields divergence- free in 
homogeneous case: (V • v='=) = 0. 

Temporary, we will not take into account dissipative 
processes. Thus, application of the canonical formalism 
becomes possible ^ which makes the analysis 
more compact. With appropriate choice for the length 
scales (~ (i+ — (Mc^/47re^?T.)^/^) and for the mass scales 
(~ M), the Lagrangian functional of the incompressible 
two-fluid model, in the absence of an external magnetic 
field, takes the form 
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Here fj, — m/M is the only dimensionless parameter re- 
maining in the system. For the electron-positron plasma 
^ = 1, for the hydrogen plasma ^ « 1/2000 <^ 1. Be- 
low we consider the latter case. The first two terms in 
the expression (^) give the kinetic energy of the ion and 
electron fluids, while the third term is the energy of the 
magnetic field created by the flows of electrically charged 
fluids. The conditions of incompressibility are assumed, 
(v^ • k) = 0. A possible derivation of this Lagrangian 
is via using the microscopic Lagrangian of a system of 
electrically charged point particles that can be written 
up to the second order on v/c, as it is explained in the 
famous book by Landau and Lifshitz [[16|: 



and x~(c, t) are incompressible Lagrangian mappings de- 
scribing the motion of points of the ion and electron 
fluids, labeled by the labels a and c. The correspond- 
ing mathematical technique is explained, for instance, in 
Refs. The equations of motion of the two-fluid 

incompressible system have the following structure: 
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dt Sv^{r] 
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where the operator in the parentheses on the r.h.s. is 
the projector onto the functional space of divergence- 
free 3D vector fields |l8). The two vector fields 
p^(r) = SCft/6v^{r) are the canonical momenta by def- 
inition. In the Fourier representation they are given by 
the expressions 
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Below, we will need the reversal relations for the veloci- 
ties through the momenta: 



{^lk^ + i)p+ + Pi, 



(fc^ + l)Pk +Pk 



_ (7) 

It is possible to reformulate the equations (Q) as equations 
for frozen- in vortices. 
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where ra(i) are the positions of the point charges Ca, 
V£i(t) = Va(i) are their velocities, 'tVabify arc the unit 
vectors in the direction between and et,. Introducing 
macroscopic averaged characteristics — concentrations 
n='=(r,t), hydrodynamical mean velocities v='=(r, i), and 
considering the special dynamic regime of slow vortical 
flows (" soft" degrees of freedom) , when deviations of the 
concentrations from the equilibrium are not excited, as 
described in Rcf.|l^, one can arrive at the expression (||). 

It is important that the variation of the action func- 
tional S = J C^dt, which is necessary for constituting 
the equations of motion, should not be performed with 
respect to the variations ^v='= (r, t) , but with respect to 
the variations (5x+ (a, t) and (5x~ (c, t) , where x+ (a, t) 



nf{r) = curl 
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where the canonical vorticity flelds are defined as the 
curls of the canonical momenta, ft^{r,t) = curlp='=(r, t), 
and also the Hamiltonian functional of the system is cal- 
culated: 



n^{n\n-}= J {(p+.v+) + (p-.v-)}dr-/:. 



(2^ 
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It is clear that in the problem under consideration 
there are two separated dimensionless scales of inverse 
length, fc+ ~ 1 and /c_ ~ 1/A, where A — ^/JI is the elec- 
tron inertial skin-depth (normalized to (i+). Since <c 
1, we may write with very good accuracy 7ip{r2+, fJ^} « 
nx{fl^,^'}, where 



Hx{n+,n 



(2^ 
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Depending on the typical spatial scale of the vortices, 
several dynamical regimes are possible in this system. 
The small and moderate wave number region, k <~ 1, 
corresponds to the Hall MHD, and in the special limit 
+ < \Q.-\, we have here the usual MHD. 

The region 1 <C /c <~ 1/A, under the extra condition 
|r2+| < |r2"|, corresponds to the EMHD [||. For the 
flows with larger typical wave numbers, k ^ 1/A, the 
magnetic effects become relatively un-significant, and the 
system ( p^ ) is broken into two weakly interacting sub- 
systems, each of them being approximately described by 

the ordinary Eulerian hydrodynamics, since G (fc) w 

1/A2fc2, G+-{k) w 1/A2fc4 < G++{k),G—{k) in this 
region. 

Axisymmetric large-scale EMHD flows. — Let us now 
consider the subset of solutions, for which the ion canon- 
ical vorticity is identically equal to zero, ft^ = 0, and 
the electron vorticity Jlj^ is concentrated in the range 
l<Cfc^l/Aof the wave numbers, where the Green's 
function G--{k) is almost flat: G--{k) « 1. Practically 
this corresponds to the condition 3 <^ k <^ 20. For 
EMHD model this is the long-scale region, where fl~ is 
proportional to the magnetic field in the leading order. 
It should be emphasized that with r2+ = the velocity 
v"*" of the ion component is not exactly zero, however, 
it is much smaller than the velocity v~ of the electron 
component, as it becomes clear from consideration of the 
Eqs. (|^) with p+ = 0. In the main approximation, the 
Hamiltonian for the electron vorticity takes the very sim- 
ple form 



nx{o,n'} 



(12) 



in accordance with the fact that the energy of the sys- 
tem is concentrated mostly in the magnetic field. The 
corresponding equation of motion is local and essentially 
coincides with Eq. (m: 



Jlj = curl [curl fl x $7 ] 



(13) 



One of remarkable properties of the equation (^^ is that 
in the case of axisymmetric flows, when 



n (r,t) uj {q,z,t)[e^^ x r], 



(14) 



where q = {x'^ + y'^)/2, we have the exactly solvable Hopf 
equation for the function uj~ {q, z, t) M: 



2lu^uj^ = 0. 



(15) 



The solution of the equation (^_5|) at i > is constructed 
from the initial function LUQ{q,z) by the shift of each 
level contour cj(^(g, z) — w along z-axis on the value 
2wt, that makes possible breaking of the profile after 



some time. Not long before the moment of the singular- 
ity formation, the equation (|l^) becomes non-applicable. 
For correction, it is sometimes sufficient to add into the 
r.h.s. of the equation ( [l^ ) the only linear dissipative term 
{e'^n/Ma)An^ , which takes into account a finite electri- 
cal conductivity cr ||^. In this case the equation for the 
function uj~{q,z,t) looks as follows: 



9 

e Ti 



(16) 



In order to justify the neglect by dispersive effects, the 
typical values of uj~ should not be too large: <~ 
e^n/2\M(j w IQe^n/Ma. With this condition the width 
of the current sheet will remain several times larger than 
the dispersive length A. Otherwise, it is necessary to 
take into account subsequent terms in the expansion of 
the Green's function G^-{k) on powers of }?k^ (we may 
neglect the term l/fc^ as previously, since fc 1): 

G__(fc) w l-X^k^ + {\^k^f + (17) 

7iA{0,fl"} w iy"j7" • (l + A2A-f...)J7"dr.(18) 

Let us consider the axisymmetric flows like (p^. It is 
useful to note that in the absence of dissipation, as follows 
from Eqs. (^), the dynamics of the functions uj^{q,z,t) 
possesses the remarkable structure: 



-{6n./SLu^)^Lu^^0, (19) 



where n^{uj+,uj~} = {l/2Tr)n^,{uj+[e^ x r], 



XV]}. 



Thus, each of the functions uj^{q, z, t) is transported by 
its own, divergence-free in (g, z)-plane, two-dimensional 
velocity field, the stream-function of which coincid- 
ing with the corresponding variational derivative of the 
Hamiltonian. The same Poisson structure governs the 
ideal hydrodynamics in Cartesian plane p3| |. 

Using the expression for the A-operator in {q, z)- 
coordinates, 

A{/(g,z)[e, X r]} = {2qf,,+Af, + /,,)[e, x r], (20) 

we easily obtain the asymptotic expansion (for simplicity, 
we write uj instead of in the two following equations) 

n^{Q,Lo}^ j uj[q + \'^{2dqq'^dq + qdl) + ...]Ludqdz 

(21) 

and the corresponding conservative equation of motion 



ojt + 2uju!z + 2A 



{2q^UJqqz + AqUJqz + qUJzzz)ujq 



+ {8qUJqq + Aujq + UJzz +2q'^LUqqq + qUJzzq)uJz =0, (22) 

where the nonlinear dispersive terms are explicitly writ- 
ten in the first order on A^ . The dissipation can be taken 
into account like in the r.h.s. of the Eq.([l^). 

In the special case when w~ is only slowly dependent 
on the radial coordinate q, but strongly depends on the 
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axial coordinate the expansion of G (fc) on the pow- 
ers of X^{kl + ky) is appropriate: 



1 A^(fcg + fcg) 



(23) 



Then in the leading order the equation of motion for 
LO~{z,t) becomes nonlocal integral-differential: 



+ 00 



(24) 

For long-scale profiles of ld this equation is approxi- 
mately reproduced by Eq. (p^), but in addition it is 
able to describe changing of the steeping regime from 
explosive Iwjlmax ~ {t* — to exponential Iwjlmax ~ 
exp C{t—t^,) after the width of the shock becomes smaller 
than A. The exponential growth of the maximum of \uj~\ 
takes place on the final stage of shock evolution (without 
dissipation) since the integral operator in Eq.(24) makes 
the transport velocity for uj~ smooth enough even for a 
very narrow shock. 

Shocks in helical EMHD flows. — Analogously, the 
helical flows can be investigated, with 

(il^)^ = ^l{x cos K z + y sin if z, y cos Kz — x suvKz, i)(25) 

that are space-periodic along z-direction with the period 
= 2Tr/K. The general solution of Eq. (|l|) for this 
case can also be obtained, since the equation of motion 
for the function Q{u,v,t) is 



nt + 2K^n{vnu - un^) = o. 



This equation follows from the Hamiltonian 



Q du dv. 



(28) 

Each level contour fl{u, v) ^ W rotates with the individ- 
ual angular velocity dO/dt = —2K^W, that is the rea- 
son for shock producing. Higher-order corrections to Eq. 
( p7| ) can be derived similarly to Eqs. (|20|-p2|). However, 
in this case it is not possible to include the dissipation 
into consideration in the framework of single- function de- 
scription (|25-26), since magnetic diffusivity destroys he- 
lical shapes of the magnetic lines. 

Shocks in Hall MHD. — If we would like to escape the 
restriction A: 3> 1, it is necessary to deal with the Hall 
MHD, the Hamiltonian of which is 



+ ^ Ji^^ + ri")(-A-i)(0+ + n-)dr. (29) 
For axisymmetric flows we have 



n 



HMHD 



{liJ^,uj } = / (oj Yqdqdz 



{ijj^ -\- uj )G{lo^ + uj )dqdz, 



(30) 



(27) where the operator G is defined as follows: 
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(31) 



The equations of motion can be written in the form 



0, (32) 
0, (33) 
G{u+ +UJ-). (34) 



Since the nonlocal operator G possesses smoothing prop- 
erties, analogously to the usual "flat" A~^-operator, the 
stream-function is smooth enough even where the func- 
tions w"*" and u}~ have infinite gradients. Therefore, the 
effect of the non-locality, generally speaking, can not 
overcome the tendency towards the breaking of the func- 
tion oj^ profile, at least with moderate typical values of 
We can suppose that with the initial data concen- 
trated in the region A: ~ 1, the breaking takes place as 



the general case in the Hall MHD model. As concerns the 
transition to the limit of usual MHD, on small fc <C 1, and 
w~ ^ '^q, \uj^ -|-a;~| <C Iw"*"!, \ijJ~\, in this case the ques- 
tion about breaking remains subtle and needs additional 
investigations. 

Estimations for the shock singularities in different 
models. — Finally, let us note that the equation ( [T^ ) 
is interesting also from a more general theoretical view- 
point. This is an example of a 3D hydrodynamic type 
system, where the singularity formation explicitly takes 
place in a finite time. However, the above described 
mechanism for the singularity formation can not be uni- 
versal for all the hydrodynamical systems. For example, 
it is known that in solutions of the Euler equation no 
singularity can form in a finite time without maximum 
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of the vorticity growing to the infinity [Q. Simultane- 
ously, the field of the vorticity direction must loose the 
smoothness at the singular point But in the cases 

considered here no of these two conditions is satisfied, 
but nevertheless, the singularity develops. We may sup- 
pose that a type of a possible singularity in some hydro- 
dynamic system depends on the behavior of the corre- 
sponding Green's function at large k. In Eulerian hydro- 
dynamics G{k) = while in the model (13) we have 
G{k) — 1. The natural question arises: If the Green's 
function has the power-like asymptotics G{k) ~ 
with a constant 7, what is the critical value 7c of the 
exponent that separates systems similar to Eulerian hy- 
drodynamics from those similar to the model (|3|) , as far 
as the axisymmetric flows are concerned? The following 
simple estimations give the answer 7c = 1. Let us con- 
sider the nonlinear transport equation like Eq.(^^ for 
the conserved quantity w(z, t), 



for s is the following: 



Kt) - -2s{t) 



ikCj{k,t)G{k)dk/2TT. (36) 



Now we note that the function —iku}{k^t), the Fourier 
transform of —ujz{z, t), is a smooth real- valued symmetric 
function with the maximum value 2 at fc = and with 
a width of distribution in fc-space of order Afc ~ s{t). 
Therefore the approximate relation is valid, 



s{t) 

[t) ~ s{t) I G{k)dk. 



(37) 



From here, in the case G{k) ^ 1/k'^ with 7 < 1, we easily 
derive the singular asymptotic behavior 



(38) 



ujtiz,t) + 2uj^{z,t) 



Co{k,t)G{k)e'''' dk /2tt = 0, (35) 



where oj{k,t) is the Fourier transform of the function 
u{z,t). For simplicity, it is convenient to deal with anti- 
symmetric solutions u}{—z,t) = —uj{z,t) that have the 
shape of a smooth step, with lim^^zpoo ^(-z, = ±1- 
It is natural to monitor the growth of the quantity 
s{t) = —ujz{0,t). The corresponding equation of motion 



while for 7 > 1 the integral in Eq.(|37|) converges at large 
values of s and therefore just exponential growth takes 
place, s ^ expC(t — t*). 
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